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Abstract: The problem of steady laminar fully 
developed flow and heat transfer in an inclined channel 
consisting of a couple-stress fluid sandwiched between 
two Newtonian fluids with permeable beds is 
investigated. The flow in the permeable beds is 
governed by Darcy's law. The closed form solutions are 
obtained for the velocity and temperature fields in the 
channel. The effect of permeability on the velocity and 
temperature is discussed. We observed that for a given 
y, the velocity increases with increasing Darcy number 
and for a given y, the temperature increases with 
increasing couple stress parameter a. 

Key words: couple-stress fluid, Newtonian fluid, inclined 
channel, permeable bed. 

I. Introduction 

The flow of non-Newtonian fluids through and past 
porous media have wider applications in many 
branches like petroleum Engineering, chemical 
engineering and such other important fields. There are 
several non Newtonian fluid models available in 
literature. Some of these models deal with theory of 
polar fluids. The theory of polar fluids was derived from 
a statistical mechanics model that assumed non central 
forces of interaction between particles. If the inter 
particle forces are not central forces in a particle 
interaction, there is an inter particle couple as well as an 
inter particle force. Due to action of this couple, the fluid 
particles will have a tendency to rotate relative to 

their neighbors. The idea of polar fluid is obtained by 
introducing a kinematic variable to model the forces that 
balance the action of the couple. 

The couple stress fluid model is very useful to 
explain the behaviour of some slurries and some 
physiological fluids. Initial studies of experimental data 
on blood suggest that some derivations from Newtonian 
behaviour of blood may be explained through couple 
stress fluid model. Based on the couple stress theory of 
Stokes (1966), Chaturani and Kaloni (1976), Chaturani 



and Upadthya (1979) have studied some theoretical 
models for blood flow through narrow tubes. 

Chaturani et.al., (1981) studied a three 
layered coquette flow model for blood flow and they 
assume that the top and bottom layers consist of 
plasma (Newtonian fluid) and the middle consist of a 
red cell suspension couple stress fluid. Umavathi and 
Malashetty (1999) studied the effects of couple stresses 
on the free convective flow in a vertical channel. Free 
convection flow of an electrically conducting couple- 
stress fluid for the radiating medium in a vertical 
channel has been studied by Umavathi 
(2000).T.K.V. Iyengar et.al, (2011) studied the pulsating 
flow of an incompressible couple stress fluid between 
permeable beds. 

Umavathi et.al., (2005) made a detailed 
study on the flow of heat transfer of a couple stress 
fluids in contact with a Newtonian fluid. Keeping in view 
the practical applications described above we have 
analyzed laminar fully developed flow and heat transfer 
in an inclined channel consisting of a couple stress fluid 
sandwiched between two Newtonian fluids bounded by 
permeable beds. The flow in the permeable beds is 
governed by Darcy's law. The closed from solutions are 
obtained for the velocity and temperature fields in the 
channel. The effect of permeability on the velocity and 
temperature is discussed. 

II. Nomenclatures 

x, y - Cartesian coordinates 

u x ,u 2 ,u 3 - Velocity in region I, region II, region III 

u b1 - Slip velocity at the upper permeable 

bed 

u _ - Slip velocity at the lower permeable 

B Lj 

bed 

- Darcy number at the upper permeable 
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Ec 
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K, 
K 2 
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T 1> T 2> T 3 

Pr 



bed 



( 



v 



J 



Darcy number at the lower permeable 



bed ^ 
Temperature at > = 2h 

Temperature at > = -h 

Angle of inclination of the channel with 
the horizontal. 
Slip parameter 

Average velocity 



Couple stress parameter 



( 



v 



JU x h 



2 \ 



Specific heat at constant pressure 

Eckert number 

Hight of the regions , I, II, III 

Ratio of thermal conductivities ^ 



h." 



Thermal conductivity of the fluid in 
regions I and III 

Thermal conductivity of the fluid in 
region II 



Ratio of viscosities — 

Pa 



Non dimensional pressure gradient , 

dp 



-Re 



a* 

Re 



pi 



sin 0 



F 



Temperatures in region I, II and III 
Prandtl number 



III Mathematical formulation of the problem: 

Consider the steady laminar flow of a couple stress 
fluids in contact with a Newtonian fluid in a channel 
bounded permeable beds as shown in fiaure 1 




Fig 1: Physical model 



The channel is inclined at an angle 6 with 
the horizontal. The regions -h < y < 0 and h < y < 2h 
are filled with a clear viscous fluid of viscosity ju x and 

thermal conductivity ^and the region 0< y <h is 
occupied by a couple stress fluid viscosity ju 2 and 

thermal conductivity^. The fluids in all regions are 

assumed to be immiscible and the transport properties 
of the fluids in all regions are assumed to be constant 
the flow in the permeable beds is described by Darcy's 
law. 

The following assumptions are made in simplifying the 
basic equations 

The flow is steady, viscous, incompressible and fully 
developed. 

The flow is in x-direction. All the physical quantities, 
except the pressure are functions of y only. 
The motion is caused by a constant pressure 



gradient 



dp 



iv. The permeability's of the lower and upper beds 



are an d k 2 ■ 



In view of the above assumptions, the equations 
of motion for velocity take the following form: 



Flow in the channel 
Basic equations 

Region I 



(D 



0 < v < 2h ) 



d u , d 



1 °P r\ 

h- — + pg, =0 



(2) 



a. 

Region II 



+ Ps , = 0 



(0<y <h) 



3 M o 3 U 



M2 ~\ 



2 



-^r+Ps, =0 



3/ 



(3) 



dp 

a7 



+ ps =0 



Region III 



(-/» < V < 0) 



a 2 



Un Up 

Oy Ox 



= 0 
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dp n 



(4) 



11 



dy 2 



= 0 



a t y 



= 0 



05) 



Similarly the equations of motion for temperature in the 
three regions take the following form: 



Region I 

d 2 T 



[h < y < 2h ) 



i Mi 
L +c — 



dy 2 



Region II 



i 



du 



dy 
\ y J 



= 0 



(5) 



(0<y<h) 



a 2 



T 



2 



dy Z k 2 



Region III: 

d 2 T 



du 



dy 

V 7 J 



= 0 



(6) 



(-h < y < 0) 



a 2 



3 +c 



3 



\ 



2 



3 



= 0 



(7) 



Flow in the porous medium 



The flow in the lower and upper beds is governed 
by the following Darcy's law. 

k . dp 

Q * = ^ ^ P 8 x 

jU Ox 

Where 

i = 1 for upper permeable bed and 
i = 2 for lower permeable bed 



Boundary Conditions 

du 1 -a 



u i = u B i> 



u, = u 2 



U 2= U 3 



U 3= U s2 > 




k 



{ u B i~Qi) 



a t v 



a t 
a t 



du 3 



dy 

du 1 d 




y —h 

v=0 



' U s2-Q 2 ) 



at v 



= 2k (8) 

(9) 
(10) 

= -/, (11) 



3 



y 



U o 

dy 



u 



2 



dy 3 



at y — h 



M 2 

<Fl 
dy 



hi 2 d u 2 du g 



dy dy 3 1 dy 



1 



a t y 



= 0 



— 0 at y — h 



(12) 



(13) 



(14) 



>T —H 




T 1~ T bV — 

T^—T 2 at y — h 
Tn — Tn at y — 0 



(^i-^o) at y = 2h 06) 



T 3 ~ T b 2> 



(17) 
08) 



H 



K 



dr x 

dfr 



K 



— K 



= K 



dy 
dr 



at y — h 



3 



dy 3 dy 

IV. IMON-DIMENSIONALIZATION OF THE FLOW QUANTITIES: 

Introduced to make the basic equations and the 
boundary conditions dimensionless: 



a t y 



• =0 



(20) 



(21) 



X — 



X 



h 



y 



y — 



* 0, 

u 



0 = 



h 

T —T 



b2 



= " 0=1,2) 

u 

_ Mi. 



m — 



* P 
P — 



1 T b\~ T b2 fa A 

In view of the above dimensionless quantities 
the basic equations and the boundary conditions for the 
velocity take the following form (Neglecting the bars). 



2 



Basic Equations 
Region I: 



, 2 

a u y 

~~^2~ = H 1 

dy 

dp cosO 



+ 



= 0 



(22) 



(23) 



>y 



F 



Region II 



4 2 2 

d u 2 a d u 2 

~2~ = 11 2 



d y 



4 



m d y 



dp cosO ~ 

— + = 0 

dy F 

where H 



(24) 



(25) 



Re 2 



2 - a 

F 



sm.&-a p 



Region III 



, 2 
d u 



d y 



2 ~ H 1 



(26) 
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dp COS0 ~ 
— + = 0 



F 



Re . 



Where h 1 = p 0 sin<9 ,p Q = Re 



o OLXA ^ ' 1 o 

F 

Boundary conditions are 

at y — 2 

d u 



a, 



F = 



U l= U B l 



(28) 



1 



dy 









-a ^ 


U 


mi-Qi 


\l Da i 






2 at 


y 


= 1 




y 


= 0 


B 2 at 




y=-l 



du 



3 



K 2 -e 2 ) 



=-1 



1 
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2 3 
m dy a dy 



a t y 
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(30) 
(31) 
(32) 

(33) 
(34) 
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= 0 



= 0 
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= 1 
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(36) 



(37) 



dy 

Similarly the basic equations and the boundary 
conditions for the temperature take the following form. 

Basic equations 
Region I: 



d 0, 



+ ec Pr 



dy 

Region II 

d 2 $ 



d u 



1 



= 0 



(38) 



K 



y dy J 



+ — ec Pr 



dy 



m 



d u 



\ 



dy 



= 0 



(39) 



where k — 



2 



Region 

d 2 0 
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+ ec Pr 



dy 



du ^ 
\ dy J 
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(40) 



Boundary conditions are 

d 9 X —H 



dy 



4 Da \ 



0 Bl at y — 2 
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(42) 
(43) 

(44) 

(45) 
(46) 



V. Solution 



Solving the equations (21) - (23) subject to the 
boundary conditions (27) - (36), we obtain the velocity 
field as 



Pr 
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2 



+ c x y +c 2 , 
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P — 



1" 1 



V 
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The slip velocity at the upper permeable bed is obtained 

as 
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u Bl = 2c 1 + c 2 + 2p (50) 
The slip velocity at the lower permeable bed is obtained 



as 



^2=- C 7 + C 8 + 



P 



(51) 



Solving the equations (38)-(40) subject to the 
boundary conditions (41)-(46) we obtain the 
temperature as 

0 1 = a 1 y 4 + a 2 y 3 + a 3 y 2 + B iy +B 2 (52) 

0 2 = a u cosh ( 2a 4 y ) + a 15 sinh ( 2a 4 y ) + a m y cosh (a A y ) 



+a 17 y sinh ( a 4 y ) + a 18 cosh ( ^ 4 v ) + a 19 sinh ( a 4 v ) 

4 3 2 

+#T 2a y +^21^ +5 3^ +5 4 (53) 

^=% 4 +% 3 +% 2 +5 5 y +5 6 -1<>' <0 (54) 
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+a A a 16 sinh ^ 4 + a A a 18 sinh ^ 4 + ^ 4 a 17 
cosh ^ 4 + a 4 a 19 cosh ^ 4 + a 16 cosh ^ 4 
+a 17 sinh ^ 4 + 4ctr 20 + 3a 21 + 2or 22 
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VI. Deductions and discussions 

From equations (1) - (4) we have calculated 
velocities U (i.e. U 1f for region I, U 2 for region II, and U 3 
for region III) as function of y, for different values of 
Darcy number Da for fixed p =-5, a =\ m-2, 6tr=Q3and 

is shown in figure 2. We observe that the velocity 
increases with the increase in y initially and then it 
decreases with the increment in y. For a given y, we 
notice that the velocity increases with increasing Darcy 
number. 

The variation of velocity U with y, are calculated 
from equations (1)-(4) for different values of couple 
stress parameter a and is shown in figure 3. We 
observe that for a given y, the velocity increases with 
increasing couple stress parameter a. 

The variations of U1, U2 and U3 with y are 
calculated from equations (1)-(4) for different values of 
viscosity ratio and are shown in figure 4. We observe 
that for a given y, the velocity decreases with the 
increasing in the viscosity ratio m. 

From equation (52)-(54) we have calculated 
temperature <9 (6^ for Region l,# 2 for Region II, and 

^ 3 for Region lll)as a function of y for different values of 
Darcy number and for fixed 

P = — 5, a = 1, m = 2, a = 0.3, K = 1, E = 1, D a =0.1 

and h = 0.6 and is shown in figure 5. We observe that 
the temperature increases with the increasing in y and 
then it decreases with the increment in y. For a given y, 
we notice that the temperature increases with an 
increase in the Darcy number. 
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The variation of temperature 0 with y are calculated 
from equations (52)-(54) for different values of couple 
stress parameter and shown in figure 6, for 
fixed P =-$ a =\m =2^=Q3^ =\e=\Da =Q\h =Q6. We 

observe that for a given y, and temperature increase 
with increasing couple stress parameter a. The variation 
of temperature <9 with y are calculated from equations 
(52)-(54) for different values of ratio of viscosities m and 
is shown in figure 7, for fixed 
p = -5, a = X in = 2,a = 0.3,K = 1,6 — \d a = 0.\h =0.6. 

We observe that for a given y, the temperature 
decreases with increasing viscosity ratio m. 
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Fig. 2 Velocity profiles for different values of Da 
with fixed values of p=-5, a=1 , m=2, ^=0.3 
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Fig. 3 Velocity profiles for different values of 'a' 
with fixed values of p=-5, a=1, Da=0.1, m=2 




Fig. 4 Velocity profiles for different values of m 
with fixed values of p=-5, a=1, a =0.3, Da=0.1 
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Fig. 5 Temperature profiles for different values ofDa with 
fixed values ofp=-5, m=2, K=1, H=0.6, a = Q.?> 
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Fig. 6 Temperature profiles for different values of 
with fixed values of p=-5, m=2,Da=0.1,K=1,H=0.6 



a' 



6 



International Journal of Innovative Technology & Creative Engineering (ISSN:2045-8711) 

Vol.2 No. 6 June 2012 



6 



1 (\ 

9 




T 

I 

y ? 


III 

/ 






8 







j- 










' i 








6 












5 








111 


4 








I 


2 


3 








II 


3 


2 








III 


4 


1 




i i i i i 


1 


l 








0 30 45 60 




90 





y 

Fig. 7 Temperature profiles for different values of m with 
fixed values of p=-5, Da=0.1, K=1, H=0.6 f a = 0.3. 
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